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TECHNICAL FIEL^^ 

uo image processing, and more 
,00011 ThepresentinvenUonrelatesin^n^^^^ 

techniqueisbroadlyre ,„^,ce projects X-wavelengtH raa ^„,e of the object, 

^ ,WCT Scanner, also called comp , .^construction from 

^iosraphic(e.g..X-raV).-g-' „^,„p,rt„{thehum^ 

instance, real-time X-ray 
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«*.rinheral vascular procedures, percutaneous 
transluminal catheter angioplasty (PTC A) proc 
procedures, as examples. 



procedures, as ca«uxih»-- 
p^duct, such as electronic devices(e.g..pnn.e<lcu».b^^^^ 

be referred to herein as 'depth layers po,w,ance, a printed circuit board (or 

evaluate thequality of theobiectCorportionthereoO^ ^^^^^^^^^^^ 
„.her Object under study)mayco,nprisev»ous^hlay^^ 

,,,.lysi.nplee.a»p.e.a*«a.id^^^^^ 

.iaesoftheboard. X*"- "^^J ;ler .he circuit board may comprise surface 

rr:;::— ;terfi..de..usres.t...^.>.e.^ 



board that may be of interest. 

■ frntn various different angles ot 

-:r;rru:::::— inradio^aphici.^^ 
ro;;:^:;:::^^---:::::— 

,„ee.,ofthe Object. Thatis.X-rays,^ybeap^.ed^^^^^ 
.spectionatvariousdifferent^S.-™^^ 
opposite the X-ray source may receive the X rays 
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geometry.X-ray beams projected ftomas,ngP^^^^^^ _^^^^^ 

asourocpoint),butonlyas.„glep.aneomeX ^^^^^^^^ 
inspection. For ins.a„ce,aooUi.a,orgnd™ayb ^^^ 

^aer U^pectionexeep. ones that passthroughahtm^^^^^ 

^„.,esonrcepo.ta.onsonep.anet.atrsuse^f«^^^^ 

pass throughtheobject and arecapturedbya^^a^ ^^^^^^^^^^^ 

CCDS). The an-beam technology was con^nonly »^ 

ofaSDbodyunderinspection. 

^.eoneWMnraging technique. conn^o^y-^^P^^^^^^^^^^^^^^^^^ 
^og^phic waging system, and the captur^ ~^ ^ 

.„HUeredBac.Pr.e^n.— ^^^^ -.cone-bean, 
seetionsoftheobjectbemgnnaged)fromt , „diographic imaging devtce m 

tomograpbyisweU— inthcart^d. — ^^^^^^^^ 

source having a cone-shaped, 3D ray outp P ^ 
i„speetio„,a2Darrayofsensorsarepost«onedon*e^s.^^ ^^^^^^^ 
,oWarrayaresynchronousiymovedre.at.vet„^e.^-Th 

or the source may be moved along a y ^.„v cases for reconstructing cross- 

u 1 ov the X-ray images are not a ID line array of 
images(e.g.,CCDU«earray)as.n*^f»b^ ^^^^^^^^^^^^ 
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h as the to-beam technology, acquired only ID projections. 
,echnology,imagingtechn„ues,suchasthefan .^he 2D intages that are captured are 



portion thereof). 

systemforstormgthelarge image data becomes a concern. 

2D image data for teconstructmg the 3D m«.ge 

■ .omoaaphyanumberotimagingtechmqucsare 
„01„ WithinX-rayabsorptton-J^^^^^^^^^^^ 
appucabletoreconstn^ctionofcross-se^on^sh^^^ ^^^^^^^ 
wasdevelopedisl^ownaslammog^h . P y ^^^^^ 

„„.ed inacoordinated fashion relative to the ob,ectWb^^^ ^^^^^_^^ 
„„.sidease.ectedfoca>p.aneleadtoab,»ed.m.^.«^^^^^^^ 
No.4.926.452,Foca.p.aneimagesa.erecon^c^^^^^^^ 

.^p,eofa,ami„ogr.hysys.m-^^^^^^^ 

i,„Honprovideson.yonecross.sec..on oftheobject^d^^ 
.ample. wheninspectingaprintedcircu.tboan.a^^^^^ 
.^tbo^andwhichcrossesaU.,n^^^^^^^^^ 

board may be desired for mspecuon of the qual y ^^^^ ^ 

,aminographic.echni,uescreateonecross-s^t.o^^^^^^^^^^ 

heightatatime, Thereafter. theboardmaybeseque^^V B^ ^ ^^^^^^^ 

.eights. Witheach.maginsi.e«tion,ap^^>>--^ 
height. For instance, in the SDXmspecaon system man 
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.■ ofrtiescintilatotcoupled With the full 
.coo*....>-nosraphyissen^y- ^^^^^^^^^^ 

oyeoonahigh-speedproducuonlm. ^^..echmoue is that *e technique 

,,,, ..couaai— ---^^^^^^^ 

.einsno...eto^ove.aao«s«.-«^^^^ 

obtamedwithmorerecenttomogrp ^tomosynthesis. Tomosynthesis 

,,,3, --ima^.tee.i.e.veU..«^^^ 

..Jhuat.nto.aminoS^^'V.-- — 

u . Asthenumberofviewsbecomeslarge geometry. A 

:rria::-o.ato.ta.e..ht.^^^ 

„tiatinsfeatuteonomosynt.e.fi.»*^^ ^^^^^^^^^^^ 
^in^mosynthesisX-rayhnagesoM^^ ^^^^^^^,,,^mes^ 

... a overlapped with difterenx sy ^sed to 

"-rCI^^^ltoss-sections. hiotHe,wo.s,one^o^-^^^^^^ 
producing a vaneiy u inspection (e.g., crob 

method, indigital tomosynthesis, tn 

combined viacomputer software. „o„osed tor combining the 

,,,, VaHousmathematica,tech»,^h.eb^^^^^ 
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■ i.r in calculation but results in different 
.a.csinmeca.ou.a.ed«oss-seouonsatd.^.e„^^^^^^^^ 
,„ard.heX-rayspot(a.o„g«>.pas-J^^X» P ^ ^^^^^ 

,„p,eBacUProjec.iontechM,»esmaybeu^« ^^^^^^^^^^^^^^^^^ 

e,^edbya«liographic«^-^2L^^^^^ 
sitnpleBackProjectiontechmquessuffersmquahty 

along the passes . . 

. «f the calculations for tomographic 

.as*eso-.aUedFil.eredBackPro,eot,cn(FBP)™e^^^^^^^^^^ ^^^^^^^^ 
.peedover.hesimp.eBa.kProiecUon.J=hmq»es. ra ^ 

,egalar^ad„gof.heX-raysourcesand.mg«^^ J ^ ^^^^^^^^^^^ 

,00. — d.peo...e..a.p...^ 
.„a.esofalobie.u„de.ins^Uo„^n>.d..ap^^^^^^^^^ 

p^ductionUnes. Becauseofoo^..peed.an<.«*e^^ ^^^^^^.^^^ 
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♦ TJor pxamole, a printed circuit 

bcardassemblylinegenemUymovesvery ^ taptoent an automated X- 

^yinspectionsysteminaproducuons^^s 
UeaUy.theinspeotionprooessismreald^u P^^^^ ^^^^^ 

^uf^turinsprocess. '""-^ "^^^ „„e minute or less. Typically, automated X- 
i„.egrttyoftensof«>ou«.dsof.o,denomt -*mo^^^^^ ^^^^^ 

.,.„speoHonsystemsseveteiyUm.t.e^^.^^^^^ 

automatedX-ray inspection systems »ese^.aUy^ ^^^^^^^^^^^ 

.peciali.ed>gh-costboa.dsthata«typ.ca..y^uta 

.peedofthemanufaetu^nsUneisnot^.-^^^^^^^^^^ 
.ystemswouldlikelybeusedmoreoftenrftheyeou 



tooughput production Unes, 

BRIEF SUMMARYOFTHEINVENnON 



process captured 2D image data (e.g.. pixeis; 

. fthenresentinvention,amethodfor 

^— :::r::;spro.ded..em.^ 

processing2Dimageda.afbrrecons,™ctmg3D™g ^^^^^^^ 
LptisesreceivingapiuraiityonDm^.^--^^^^^^^^^^^ 

*"^-'^^^^°'"^l"r^.e.useaconiugategradientaigori.im.for 
herein, such iterative processmg may, for examp 

^nstructingtheSDimagedata. ,„,«.,p,esent invention, a system 

100.0. - 

,„rprocessi„g-i»agedata.rr^r.tn.m.-^^^^^^ 

..stemcomprisesameansforrece-vmgapU^.^ ^ ^^^_^^^3„^,, 

view, and ameans for iteratively processmg the reccved 
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data therefrom. Again, u 
I^esotvi.w.Thesysten,f^«con^ns^-^^ 

— '"^"^ ,^.^aa.v«.efea..e..a.ec^ca: 
,0022, Thefo^goinBhasoutUned^^^^^^^^^^^^ 

described hereinafter which fonn me subject otth^^ _^^^^^ 
::eci.edbvthoses.«edin*ear.d,a.,^»^^ 
^,,ereadi.y»tiU.edasabasis.r.o.f^.«^^^^^^ 
.^epu^sesofthepresenti— ^^^^^^ 
s„ehec,uivale„.cons™ctionsdo„o d^ft- P^^^^^ 
i„tt,eappendedc,aims. Thcnoveifeet^*^*^ 
i„ve„.on,bo.has.oi.sorga„iza.on-d»^^c P^^^^^^^^^^^ 

„eswiUbebe«erunderstood^.>^^^^»J^^^^^^,^^^^^^ 

e„.„ectionwi*«.eaccornp-y™.«i";j-^^^l,,^,,,Uononly»dis™^ 
of the figures is P«.vided for mepurposeofrtostra 

:!di"ionofthe«tsofthepresehtinve„«o.. 

BRIEFDESCRIPTIONOFTHEDRAWINGS 

f .V,. r^resent invention, reference is 
,0,3, .ora^oreco^pUteunderst^d^go^^^^^ 
„„w«adeto«tefo..wi„.Oescrip.io„s.a.en.nco.»ct,<» 

which: pvamnle 3D tomography 
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. f„, constructing 3D images in accordance with 
ZDradiographicimagcsthatmaybeprocessedforr^onstru 

embodimcntsoftheprcsentinvention; 

w 1 „h,ect that has been divided into 125 (5x5xb) 
[00251 FIGURE2Asho«sacab.calob,ectthath 

identically sized cubic voxels; „,p,™E2A- 
,00., P10UKB2Bshows.«olaye.ofvoxeU.rtheob.ecto..™ZA. 

mnle tadiographic imaging system that may be 
10027, FIGU^^^^-'^^Ttrilmaybeprocessedfor 

,0023, nOUKBS3A-3I>sho.examp,ecross.s.tiona.imagesth.maybe 
p„ducedbya3Drecons.ructionprocessforasampleob,eet.. 

[0.., HOX:KB.sho.sanexamp.eprintedcireuitbo.dass^b.ytha.maybe 
>^ec.edinaccordancewi.hembodim.tsonhepresentinvent.on., 

system; 

^a.^assignedto.hegeome.ricalprmnt.vcsof.heexamp. 

,«,3„ HGmB7sh„ws2Ds.icescombi„edintoatilethatrepresentsa3D 
volumeofthe example test sceneofFIGURES5A-5D; 

embodiments of the present invention; 

conjugate Gradient (CO) algorithm after each step of tteraUo 
scene otFIGURES5A-5D; 
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^nle «»h illustrating the observed migr^ion of 
1,035, FIGURE 10 shows an example ^'^^ ^..^thm is »sed for 

«..solutioni«to„on.phvsicaUnesative)c.ss.^^^^^^^^^^^ 
proeessingtheexampletestseeneofHOURBSSA-SB, 

. ,„™hthaliltetratesanexamplel^ttrvethat.s 

10036, FIGURB n ^o""^^".*"'' ^.cGalgorithmisusedtorprocessing 
3ssooia,edwi.h,he3Dto»ographioreconstruct,on.hentheCG 

theexampletestsceneofFIGURESSA-SD; 
ttaditional tomosynthesis technique for 3D tomograp 

of FIGURES 5A-5D; „ ,Uces of FIGURE 12A 

.0038, PIGURE12Bshowsati,e.hatccnipnses.he32shcesof 

^tHlogbrightnessnormalizationhavingbeenappUeathereto; 

rr^CrseeneofFIGURES... 

,0,0, HGURE13Bshowsa.ilethatcomp,isesthe3.sheesofHGUKBI3A 
...egbriUessno— onha^ngbeenappheathereto; 

„.n.e hidi-level block diagram for a system 

,041, ™"-^^^"::":^tfItwim entofthepresentinventio. 
toplementingaSDimager^nstrucnonprocessot 

„„le «aph of points in 3D space for which a 
,042, H°^>'*°'"'":^Xlp„intlmaybecomputed. 
Laplacianoperatoractingonthe^pomtsando^iteredat 

I^ETAILEDDESCRimOKOFTHEmVE™ 

- ..ewoft.abo.^0...;— 
development -<ia.otofappUcationsmmany.e.^^^^^^^^^^ 
nondestructive inspection, airport securuy. e^^ 

technique thMmaintainshigh,uality3I^re^—n ^ ^^^^^^^^^^^ 

«.„ghputof.he3Drecor^-on-^^ 

each Other because achieving high quahty re 
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. f^m^ sisnificant computational processing of the 2D 
,„,nnnt»«rofZDX-.V-Sesan^^^^^^^ 

— ..,,«,e.fo.3Dto.o.aphv.tHe30X 

,0., ------ ;:reolon..nsea....ee,ect.onic 

,.peotionsystentn>annfact.e.byA.»^T«^^ 

industry for the fastx-ravinspecuonofpn^d..^^^^^^^^^^^^^^^^^^^^^^ 

,^es.T.e5DXinspections^te.,snoc^jeo^^ 

hecauseoftherequiredhighthrcughputctthesy^^C g 5OX system operates utiUzing an 

production lines onwWchitis.ypica«y.»^l—^-^^ 
Log.Wtec.u.,ueinwMc.o,.tvo^e-^^^^^ 

laminographicmemodbasedonflxedp l^an^^--^^ ^^^^^^^^^ 
.suit, reeo^tn^ctedimagesofvoiume layers ^ff^^-^^^^^^^^^^^^^ 

„.er,apping.ayers.dHavelo.re..tto^»---^^^^^^^^^^ 
p,medcircuit.oardsund.r^o^--^^_^^^^^^ 

— — *""^"'"^,......se.stforeapt.n...age 

--^*"^":;rrI~<iatatoreeonstruct3Bima. 
datafbranobieetunderinspecuonandprocesstng*^^^^^^^^ ^^^^^^ 

datatheref.m.Forexamp.e,araaio^P~^^^^^^ 
e„ne-beamtecl.i,ueforeaptuHt.g2D.t»^^^ - ^^^^ 
.^^gedatamaybeprocessedusmgama^— ^^^^^^^ 

in^age data therefrom. However, theFBPt^taq ^^^^^^^^^^^^^^^ 

.^agedatainsuchSOtomographtcsy.^ -^^^ 

orderofaeqniri«gthe2D.magedata. F^'-'^l^ , .^e various different angles of view) 
directions from whichthe2Dimagedata,sa«,uu^0x_ ^^^^^^^ 

.everyorderlyassigned. Morespeeifiea..y,aconsta^^ e^^^^ 
3e,uiri„gthe2Dimagedataisgenerallyre<,u.red.Tltereare 
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. , vervregular geometry for implementing the 
i„.plementingFBP,butaUofthemethodsre.uireaveryregu 

X-ray source and sensor array. 

obiectoanbe acquired atvanous f "f^^, ^ ^,„„^ pBP techn^ue for 
^..eceiveme 2D image dat^andmay ^^^^^ ^^^^^ 

^„s«cting3Dimageda.fto.the2D,n>a^^-2^ 
ge„erallyre<.uiretha,*eobjec.u„der.nspeet,o„ e™^^^^^^ ^^^^ 

:,„..epso.— ...^^-.^^^^^^ 

aifteren. angles of view. For ms^ce, a FB ^ 
^.gingat 500 different a„g.esofv.ew*at»^^y^^ 
inspeo.ion»aybere,uiredrcbe»..c^»^^^^^^ 
i„agedauattt«.po.i«on)>e„berota.edbyo„ de^ 

.„ased.»d.enber„.edbyon^»o.degre--^ 

.on. Bach««e*eob.e..K>«.^.-^ eh confirmation, me distance tbat^e 

^lativedistanee between., and theX-«y . ^^^^^^^^^^^ 
object is translated a,ongtheZ-ax,s,sW^^y^^^^^ 
translation. Accordingly, the acqmsmon of 2D unage 
tobeperformedaccordingtoaveryregularorder, 

. .w .he 2D images be acquired in such 
,0047, Havingsnchastringentrequ.^^^^-^^^.^^^^^,^ 

..^arlyorderedwayis— tage.s«^^ 
.ppUcatlonsitn^ybedesirabetoacqune^^^^^^^ 
ofview that are notnecessonly regularly ordeMan^P^^^^^^^ 

reconstructingSDirnagedata. For instance, ..maybe^-^^^^^ 
inspection toafirstposition for imaging(.«qu.s.Uo„of2D^^^^^ 
^ew),thenrotated.eobjec.byonedegreefornna^g«-w 

..o.ectby.o.her20de.e^^^ 
desired that the acqmsition of the lu imdg 
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^Vuc.echm,uestha.r«,ui«2Dimag.datatobeao<,u,«d 
ecooMi„gtoa«ga.»oMer,s»chasFBP.^^2,, ^^^^^^ 

,„„Ubo=^isge„.a>>vave.plana.^--^^^^^ 
,„^.bca„.i.much neater *a„. *— ^^^^^,„,,,,,,„«,o^ 
ei^uUboardare— lysogreatt>..X,a^^^ ^^^^^^^^^^^^ 
,e„g«>a.d«id*isft»itless. That ..X-ray be^J^^ 

,o.d.Thus.meuseM2Din.geaa.a,s»^»«M^^^^^ 

,..edoi.u.boaM.Tb..i™as»goa^^^^ 
i^,edaUsui.ab,.torprocess.n s>sno»^^ ^^^^^^^^.^^ 

„».be.of2D images a.ai«e.eata„S.esot~ 
^onsm.eeng.be3Din>ageda..ms^^-"^^^ 

....U.g— ed30»ag _^,,.„...ebas 
J0049, Ano^erdisadv^Ugeof— ^3 

f„,p— g. Forins.ance.a«admo«^^^^ 
J„^ctioa»ay.,u..abou.500d,ffe.a^^^^^^ 

differentang^eofview. ^ '"^^ .Rations. That is, .o«ge of su.h a 

^.™„be,onOi»ages.ay — ^^^^^^^^^^^^^ 
p^eessingsuchalarge number of2D.mage«my 

> ^,„ he acceptable for many apphcations. 
^e.colong)toheaccep 3D tomographic technique ma. 

iOOSO, -ewoftheabove.ad----laregLorder. further, a 
aoesnot require ac,uisitio„of2D.mageda.mJ^^^^^ ^^^^ 
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tot allows 2D images ,0 be acqoued ftom aA.«^ .^^^^^.^^ ^^^^^^^^^^ 

Thatis.eettamembodmientsoftheprese„t regular o«ler. Certain 

C.,,X..y)i™agesl.eac<,ui.eainaeco..«w^^^^^^ 

.„.oai™e„tsoftl>eptese„tinve„.on™^^^^^^^^^^^ 

of acquiring the different angles of view. 

,0. — e.— 
SDtomo^aphicmethods that are able to .gntfieanUy^-^^^ 
applications, suchasthe above— nedtomograp^^n^^^^^^^ 

^ati^beeauseaBDlontograpbictecluu^uetsp— 

„.am>er that minimi.es the amount of such 2D g^^^^ ^^.^_^^^^ 
^.en application. For ins.ance.inag.vena^U~o^^P^^^^^^^^ 
.03^.m„cb.weranglesofview.hanare«^o^^^^^ 

^°'Ta:r:r:c:ii.nspectionce..^ 
rrri;=:^=— 

,„0S3, Bmbodimentsoftbepresentinvent.onma,ben^^^ 
^onstruotionofSDimagedataforabroadclassofseom^c.^^^ 
so^ceandX-raysensorstbatmaybeimplementeamat^ ^aph.0 ^^^^^^^^^^^^ 
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for example, sevenl times compared ^,3 of 3D U-mographrc 

Embodiments otmepresentinvenuonmayb^n-i ^^^^^^^^^^^^ 
systems. Tbatis,em— of mepres».m^^^^^ 
a...ore.„s«u.3.tmaseda« — — ^^^^^ 
„.,beemplo,edfo»aptun„.*^Dm>a^^^^^^^^^^^^ 
parallelbeam or cone-beam geomeme^ N^^ ^^^^^^^^^ 
fl.e2Dimagedata(andUscorresp<md.ngh.dw^^^ embodiments, a simple 

utilized with a specific geometry of a sys _ ' algorithm enables the 

forther below, 2D radiographic images v .8-, 

Lohic images are procebbcu 
eapturedby aradiographic imaging system and^c ^^^^^^^^^ 

reeonstrnc.3Bimagedata(e.g.,cross.s.«on^^^^ 

^ograpMcimagingsystem suchasa. -b^^^^^ 

eapw« the2D radiographic m.ageso(anob3^t^ _^^^.^^^ 

r„hicimagi„gsys.emoutputsthe2D^aa.^^^^^^^ 

en.bodimentsofthep.esentinvenuo„>ec^^^^ 

toputtoanimagereconstntctionprocesso detail below, such image 

datafortheobject under inspection. /.de^n^^^^^^ 

reconstrnctionprocessorimplemen.— 
gradient algorithm, forprocessmgthe 2D miagest 
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simitevaluesotthesolderjomt.Thevan . inspector) to detenmne 

«,enbeevaluatedbytheauton.atedinspecuonsys.em(orbyah 

whether the solder joint is of acceptable quahty. 

.he thickness of solder material (which is typically a 
10057, ^-"-'''^•JltXanautontatedinspectionsyste.*^^^^ 
eombinationofleadandtm)maybemspectedby 

a„MysisofX-rayimage(s)ofthesoldermatenal^^^^^^^^ 

„a.erialforn,ingtheX-ray.mag.W- J^^^^^^^^^^^^ 
eo^spondingto lifter shadesofgray» ^^^^^^^j^.,.„^ 
O,laok)asthethio^essofthesolderrna^al„^^^^ 

.olderwillhaveagray level that .slessthantheg^yl^^ ^^^ ^^ 

.Ider. Theimageofthethinsectie„w,,lapp^.^^eah^^^^^^^^^^^ 

.Uethictesection. ™s convention ts typ.cally ^edur^ ^^^^ ^ 

------°-t::r:::r:v;-.ofthein>ageofathic.r 

section of solder has a gray level ^^.'^^^.^^^^^^.^fonowedin fihn radiography 
sectionofsolder. '^^ — -^^T;! ^iL convention ntay be implemented 
where the X-ray images are recorded on X-ray BUn 
with embodiments of the present invenfon. 

aigitauomography system. Var,ous3Dd.g,tal«.m^^^^^^^ 
.,eart,anyofwhichmaybensedm<...»^on^- 

thus example imaging systems are only D y i„«ge data for an object 

..ayf^mthcinventivesystemandm^ ™— 
^derinspection-ExampleSPtomo^physys^^ ^^^^^^^^^ 
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inspection ftom the captured 2D image dat^ 

™, FIGURE 1 shows a schematic representation of an example 3D drgltal 
,0,591 ^^'^'r .^^^«,entbodimen.sof.heprese„.mven«on 
.omographysystem geometry whtch maybe us^^^^^ 

M„respecificaUy,FIGmElshows»cx».p^-e^_^ ^^^^ 

_.0emrtsX.ys.3.a.ne«^^^^^^ 

acircuithoard)tobem.aged^ ^eXra^^^^^ ^^^^^^^^^^^ 

„nltiple views, X-ray source 20 may travel m ^^^^^^^ 
traveling inacoordinatedcircularpath. Note*atX-«ys^2^^^^^^^^ 
^^esurraceof^X-raytubewithat-electronb^^^B y^^^^^^ 

beam.theX-raysourceise«ective.ymove.^^*^^ 

„„.actua..ymove(bu.mayinsteadbeposrt.onailyfix.^ ^^^^^^^^^ 
movement maybe achieved throng other techm,ue,^H.^V 
(e ^ to achieve imaging of the object 10 at vanous dtfferent ang 

Whileanex.plec.rcularpa.hofsot.e20..^.^ 

nOU.BU.shou.dbeu„derstoodthatm.^o.er— ^^^^^ 

nrulUpleradiographicimagesthatmaybeusedf^recons^^ 

„.jectlO,andembodime„tsof.hepreser,tmv^uo„.^^^^^^^^ 

,,e„era,foracomp.ete3Dr.— 

object. ' i^^objectlO. However, it provides 

™.hematicaUyoomplete,3Drec„nstrucnonof.^^^^^^ 

sufficient data for practical applicatton grven the ^^P* J , „„enient 

ofin.eresttobereconstruc.ed.lnaddition,*escan^«^ 
forscamri„glarge,nearlyplanar. objects. Large planar objects (eg , yp 
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!;:ig«,e.u.een.sensonoben,ovedco»p.ete.y.o™d.ts...vo,™e. 

,0061, Ofoourse,variouso*eroonfigurationsofadigi,.ltadiographici»agmg 
100611 Of^"" . 2D,adiographicimag=sofanobjec.l0maybe,mplemented 
system operable «>oapWred,g,.a. 2D rad.^^^^^^^^^^^^^^^^^^^,^^^^ 
»ifl>embodime„tsofthepresentmvenuon,melutogw,th ,,,bastha.ot 

u ftl7R223 While a cone-beam radiographic imaging system, su 
Patem Number 6,178,223. wmie . ^ ^b„diment of the present invention, it 

FIGUMUs u^ed in conjunction withapreferredembod^men P ^^^ 

• hat the oresent invention is not limited to a specihcconngui 
shouldberecogmzedthatthepre ^ ^.configuration of a digital radiographic imaging 

described hereafter for efficiently reconstructmg, from the 2D radiogrp 

data (e.g., cross-sectional images) of object 10. 

10,6Z, ConventionanDdigitaltomo^aphyapproximates^objectbydiscre^ 

„ . "Voxels" (or "volume pixels") are well known m the art of 

volumeelements,calledvoxels. ^"^'^ J i„ general, a voxel is flie smaUes. 

toage processmg, and are commonly usedm3D.mag^ng. g^^^^^^^^ ^ ^^^^^^ 

disti„guishablebox-shapedpartofathree-dimens.na^.magJ— 
.aingdepth.o.imageusingasetofcross-s.n^.»a^ 

Thesecross-sectionalimagesCorshceslaremadeupofpucds^^P^^^^^ ^^^^^ 
,.els.onesUceisreferredto.i.^~^^^ 

And, the distance between any two slices IS reterr stacked in computer 

.al-worlddepth.Thevolumetricda.asetispn.esse^whensl^^^^^^^^^ 
^emorybasedoninterpixelandintersUcedistar^^^toac^^^^^^ 
volume. Next, additionalslicesmaybecreated — ^ 

^^^^ 

10063, Forbetterunderstandingofee,.ainprincipalscommonlyutili.edm3D 
.^ographyandwhichmaybeutili^dincertain— ntsofthepresentinvention 
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,*eexan>p.eradicgrapV.cimagingsys.e»c^^^a^ ^^^^^^^^^ 
sensotssuchasa^ySOinnGOTEl. 

.„m„lete 3D image is not necessary. For example, 
,00651 Formanyapp~.aeomple.e3D g ^^^^.^^^^ 

,.^«Uono.aaouMe-siaeacircnii.oa.(.c..^^^^^^^ 

,anesor— .-^c.^sj^^^^^ 
detennine solder jomt quality. FIGURt^ 



, • ,„FTOURE2C,anexampleradiographieimag.ngsystem200, 
,0066, — ;^°7,'^.„a,178,223,,ha.mayl,eusedtor 
.ehasthatdeseribeamoreMlymU.^-^' ^^^^^ 
c^turi„g2Dradiograp,nc.magesofob,^0 ^^^^ 

^liographio imaging system, souree 20 mcludes ^ 
airec.ingapluralityo«-ray.eams(w.c m^.^^^^^^^ 

.scribed — too^ectlOatva^^^^^^^^^ 
,sedintheex»plesystem200^F^^ 

differentanglesforcapturmgs»ffic.ent2Dm.g ^^^^^^^^^^ 

„.iec.lO,itsHou,dberecog„i.edtha.va^o.o*«^^^^^^^ 

for simultaneously or serially generatmg X-ray beams tha 
:;Ltanglesorviewmaybeusedino.bercon5guraUons. 
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1. .firstX-ray beam subpath 201 and second X-ray 

bean, subpa.h 205 are ones of m»yX-raybe^sp ^^^^^^^^^^^^^^^^,^^ 
e„m»atorgrid2U. TheremaipingX-raybea^^a^ ^^^^^^^^^^^ 

X..ybea™subpa«.205pass«^ugbob.^l^^ds«^^^^^^^^^^ 
„,„*detec».anay30..twmbe..^o^*at ^^^^ 
X-raystravelingal^gfirstX-raybea^supathO ^^^^ 
„itHnobieetlO;.«.«mep.bof«>efi.tX-»yb^^^^^^^^^^^^^^^^^^^ 

.eco„asUce,235,an..b.dsUce240.P— X™^^^ 

s„bpa*201may— te«>reco„s— oa -^^^^^ 
3ecUo„a,i™ageco„esponaingtodep*Uye.2^tor^n^ 

^ rav beam subpath 205, the 2D information 
10069, WithrespecttosecondX-raybe^ J 
p.videaby™227may — ^no a ^^^^^^^^^^^^ 

L..rueteacross-sectio„.imageco— ^^^^^^^^ 
voxel465inareco„s«uctedcross-sect,onal.^ ^^^^ 
instruction of a sixth voxel 270 m a reoonstnKted cross 



depth layer 240. ,u too are 

K „h 975 and fourth X-ray beam subpath 280 are 
,0070, AthirdX-raybeamsubpath275andfo ^^^^gX-ray 

.oofm^U-^.— ""^^ " 
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. *,r. 285 are not shown for the sake of clarity and 
.ean-suhpathsenu-^ting^mseoondapert. 28^^^^^^^ 

explanation. So™eot.heX-raysthat.rave.a,on X„^^^^ 

.hpath^SOpassthrou^obiectlOandstt^ede^^ 9^a^^^^^ ^^^^^^^ 

,„.ew.ths„hpaths.01ana.03,thel«— -^^^ 

«velh,ga.ongthndX-ravbean>subpath275doesn^^^^^ ^^^^^^^ 

ehiectlO;rathertheintensi.y— .ontsanag^ega^^^^^^^^ 

„ts all plane/sUces between coU»««>r^ ^^ ^^,^^^^^^ 

inoludingtheplane/sUeescontainmgvo^^ ^ 

de.eo,or291 by X-ravstraveUng along fourthXr^y^^^^ ^^^^^^^^^^ 

.,slnglepolntwlth,nob.0tl0.ath^*e«^^^ 

output byradiographic imag.ng system 200 to an . g^^^ ^^^^^^ ^^^^^^^^^^ 

aeseribedn,o.«.lybelow^.-— — ^ 
p„theteceivedp«eldaubyf^^^^^^^^ 

voxel from all of the detectors that detect X r y ^ 
eo.npletelyorpartlallycolncidentwiththatpa.c^-^^^^^ 
...epnrposeof reconstruction -exam K^^^^^^^^^ 

„y beam subpath205 inaflrstpuel output the^by a ^^.^^^^^^^ 



reconstructing the 3D voxel data. 



reconstrucxmg uiv- -^-^ 
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applicable to any type of system wluch ac,««es 2D 

,0073, FlGUKES3A-3Dsho«exa»p.e<-ss-secUo„alin.agesfl.a.n,aybe 
[007J1 no ,«„u;^PtlO For illustrative purposes, object lu 

p^oducedbyarecons— process forasample obiectlO^Fo^^ ^ ^^^^ 

bowntaHGX;i^3AHas.es.pa««nsin*eshapeofa^^»^S 
e™beddedwi*infl.eobject 10 in tee different planes (or dep* layers) 



respectively. 

„07S, Similarly.FIGURE3Cshowsasan.plecross-sectio^Hmage^' 

high-quality image of the cross-secuon <u> 
below. 

„07«, FlGURE3Dshowsa««.plecross-sec.ionaUmage(or".omograph")of 
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■ «120 Again preferablyblurredregicnlZZand/orttseffectsonthecross 

taagel20ismmmuzedmtherccons™ f.heoresent invention below, 
e^ss-sectionas described wifl.=mbod.mentsot.hepresentm 

,00,,, W..eHOUKBS3.-30s.ow.e.an.,eo,.UO^^^^^ 
.O..Ob,.d.Oe«..n.V.e— ed.— ™ -'^^^^^^^^^ 
a^tobjectlOisintendedasageneneexa^leto Uu^«^^ 

a^«..ye.ot=.obiec.ascross.sect.on.,n«^^I^^^^^ 

having any„umberofdep*.ayersn.aybe.,«pec.ed^a^^»y^^^^^^^^^^ 

,e reconstructed as cn-ss-s^tionalimage.. For exan,p.e.a^^ 

..pUn.en.ed.i.b.a.ontatedi„specUonsys.».^^^ 

circuitboard. FIGURE4depictsa«examplepnn,edcn^^ ^^^^^^^^^^ 

circuit board assembly. 

,o™. Tbus.HOUKB.provide.aconcretee.anr,eo.an^^-»^^^ 

,.,..edin^cewitbcn,W.e.so^.^^^^^ 
thepresent invention .enotintendedtobei-mt^isoe^^^^^^^^^ 
.„.,eso.ac.uitboard.b„.n.ay^.e.^^^^^^ 

^econst^ction of 3D images from cap««ed 2D 'ad ^^^^^^^^ enviromnent 

^itboutlimitaUoninspectionofvariousothertypesofpr^u^sr^ 

..,„aa.ycontro,useinantomated.nsp— 

eon.ainedtherein(e.g..sec«r,tysy^emsfortn^^t,ngP Various other applications of 

the image reconstruction process oi cu 
those of ordinary skill in the art. 

,00., .™^.----^''r"rdr::— 

.^ponentssolderedtoeacltoftwo — su^^C^^^^ 

ontoafirst surface, component 408 tsmountedontos^^d^ ^^^^^^ 
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embodiments of the invention. 

surftce mount device. One plane (or depft lay^). P ^ 
s.faoeofci.ouitboa.d400passes«*e.ea^andso|d.^^^ 

and 406. Another plane (or depth layer), depietea y 406 and 408. It 

of.reuith„ard400passes.hro„^d>eleadsand»l^«^^^^^^^^ 

.o.d.ennder..<»d.hat-etw„e.a.pe^^^ 
e.a»p.e,inotherexa»plesanyn«»Wd^th ^^^^^^^^^^ 

.^er,whUeva.oustyp.o-^^^^^^^^ 

inthisexampleasbemg.ncludedwrttan , be understood that m 

e„„po„e„.s402,404,a„d406falU„gw..h,ndep.h^y^ 1 ^^^^^^^ ^^^^^^^^ 
..eintple.ent,tlonsditrerenttype.ofsolder,o.^^— 
.ayersotthecireuitboardtwhioh^aytncreasethen^n^o e^^^^^ 

^pecUonsystent,On.yonedepthlayereorre.o^^^^^^^ 

..owinthis example for sinrphcityandeaseof — ^ 

„,,,.,ayer..naybeo,.«^t— ^ 

certain appUeations to have a 3D rmageofe ^^^^^^ 
provide sufficient infornuttion.odetennme.hat each P^^^ ^^^^^^^ 
„sefulinfo.ma.ionforinspectionofthevanoussolder,omts(e.g. 



such solder joints). 
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•♦v,CTr,TTRFS3A-3D. Once produced, 
3D images of ^.solderjomuofboard 400 oro P^^^^^^^^^^^^ 



to the user. 

be acquired, stored, and analyzed by the compu • ^ ,,..^130 encompass computer data 

i^ge. For instance, i„cer.ainembod,me„.smer— 
^,es)™a,.eais,a.ea.aus«foH„sp.^^^^^^^^^ 

computer for quality control. 

conjugate gradient method (CG) .saw ^^^^^^^ „ , ,3 

,argesysten«oflineare<,uat.ons.CG.sefl«>« J^^^^^^^^^^^^^ 
,„^ow„vector.Msatoown vector. »d.J^^^^^ 
<.„on-negative)ma.H.. Iterativeme.0^^^^^^ 

conjugate gradient method.sdescr.b«Jmorefi.m^m y^^ ^^^^^^^^ 
article titled-Anintroductionto the conjugate Gre^en l^ 

Pi.tsb«rgh,PA15213(August4 1994 * ^^^^^^^^^^ ^^^^^ 
.ferenceCandacopyofwhtchtsa^^^^^^^^^^^^^^ 
application at htffi;teffl=SLan!t^ai^ 
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hereby incorporatedhereinby reference. Tl«..s.Etoa^l^^^^^ ^^^^^^ 
proposed using the CGmeftod for proeessmgcap^redlDm-aged 

image data therefrom. 

,0,85, However,itwasno.read«yapparentfton,*ewor.byPicoo.ominiand 
. CoL«rodwasexte„dab.e,o3Bin,agereeo„struction. Tbe re^.an.a.on 

irrZaLnbe.ow,e„ab,es.eCGn,e*od.^ 
term Y Ml ^ t,<luaiion formula but 

^Menr. I- we„ ^own ^aUbe regu.ari.aUon, even ^r.^ 
^Uedforspaceswi^differentdi— w^^^^^^^^^^ 
.ownintbeabove^icieofP— -^^^^^^ 

after some number of uerauons (before* CGm^ ^^^^^^^^^ 
resuUingSDimageswasonparwimares.andardm-'^^^^^^ 

SOin-agere— c^n^P—^^^ 

take for the 3D ease (to be practical) an applicable for 3D image 

AS discussed fiHtherbe,„w.Ihavediso„vered«,a.*eCGm^*^^^^^^ 

_nstn.ctionfrom2Dimagedata«i.hanacc^t^^u.b«^^^^^^ 
,„^ity 3D images. More specifically,asdes<.bed^^b^w.-^^^ 

coefncient....ao<^tab.enurnberofterat.<^^^^^^^^ 

,mye.periments)provid.ave..o«^^^ 

n„mberofiterati„ns(e.. 

computational simphcty and the 0-own) r ^ ^ 

superior to other tomographic techmques for such^^^ W^^ ^^^^^ 

printedcircuitboards,forexamp,e,thatre,u,rev^jeMy^P^^^^^^^^ 

fl»ughput.Fur.her,theCGme.hod,wh»app..ed»3Dm^^^^^^^^^^^^ 

3S good asit does forthe 2D case(descrrbedmtheabove^^^b^^^^^ ^_^_^^^^ 

Ar^d it is onpar with the standardSDimageprocessmg methods. UkeF 
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for processing). 

,00861 TheCGn.efl,odbelo„gs,omefamnyofgradien.-basedalgorithn« in 
10086] ,^^,„ti,™ a straightforward mathematical approach to the 

ti,™ the X-rav image can be written as 
voxels in a vector j:, then the X-ray miag 

,,nXrav image (log value of pixels of the X-ray image are 
where/is a vector representmg a 2D X-ray mtage ( g dimensionality of the 

organized into the vector^andPisamatnxofpro^ect-ht^^^^^^^ 
vec.orx.„«isn.uchlarger.ha„thedimenstona.nof«,e™ 

mstead of Equation 1 above, one works with a set of equations 



PX = J, I' = / , 

oalculationbelow. Then, .hes„lutionofEq«atio„2ca„be expressed. 



.=M (3) 



for L2 metric as: 
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lOOST, The„o„eoandifferentiateE,ua.io„4by.andget*cfoUo«ingUnea. 

equation: 

|p,'(p,x-/,)=o 

' ,0„8», Or,re^«ngE<^on5in*,*e^d».fon„withtheu«,=nownonfl.e 

left side: 

r . (6) 

10089] Finally.Equationecanbewrittenasastandardlinearequ^^^^^^ 

(7) 

Ax = b 

where 

AJp^. (8) 

1=1 

This means that matrix ^compnses(lOOM) elements ,o* Qbvtes of memory, 

lines. 

„„„, TheCGmethodprovidesasolmion. This iteraUve metbod needs only to 
, ,« „f^-yp'Px,whichca„beconducteddirectlyftomthephysical 
perfonn the calculation ofAx-2^^f, • 

p In certain embodiments, 
formulationottheproblem,withoutevenstor.ng.hematncesP,.P.Inc 

lomiuiat r .t,«fi„/'Je durine run time of the image 

theP^ termintheaboveequationscanbecalculatedonthefly(i.e.,dunng 
theP termininedu ^ ^:^^tQ the P^ term is a mathematical 

processing). More specifically, in certam embodiments the P term 
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^„se„.ationof.hebac.p™ieo.io„. and thus c^beea^lyceloulatedinacco^^^ 
techmques known to those of ordinary skill in the art. 

ll,odisprete^becauseitU.heonlyiteraHven>c.h«lo»re„tlykno™(n..heory)*at 
guarantees the convergence foralimitcdnwnbcrofiterations,„. 

100,3, F„re.arnp,e>esteepes,descen.(gradient»ethod)nughtbeused»^ 
„lblentsbecausetbesradlcntinB,„aU„ns3and4.oveeanbec.cu,at«^ 

explicitly: 

slowconvergence.comparedw.ththeCGmethod. 1 ,<,^ided by 

to the CG method. 

algorithms, when the COnretbodprogressesalongnewdireCons of descent .tusesthe 
knowledge for all subspaces that participated in thepreviousiterauons. 
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cross sections, as examples, inemaui • .^.n methods are widely used in 

-rrT— ^^^^ 

introduceda/jrioryknowledgeanau v«nwledee Among the standard 

curve represents *e best value for the repolarization parameter, y . 

,0(«6, T1.est„poriteriafor.heCGi.erationsn.ayno.bederiveddireot,yftomme 
100961 ines p .^.^t^e exact solutionfor the particular. magmg 

„hservedrecons.rucHo„oteross-sect.o^^b«^^^^^ 

s,ste.heingusedn.a.be««Wo^>^^^^ 
the "model" X-ray images obtamed ftom the reconstruc 
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,essman5foroth«i»agingsys.ems/.as^,a„da„y-hn«mberofCGiten>«o„sisin.^^^ 

be within the scope of the present invention. 

,00971 Thereisavarietyofpossibleregularizationstotheaboveproblem, which 
enable the gradient algorithn, to succeed in reconstruction of 3D images 6om2D images. To 
select a regularization, one should also follow restrictions imposed by the nature of the 
„animi.a:na,gorithm.^eCGalgori.hm.initsst.dardform.assumes^top^^^ 
U„e..nom,egative symmetric matrix. Asaresult.thebestwebave„hand^r*^^^^^ 
is the Laplacianoperator, which suppresses thehighftequencycomponentofthes^tauo. 
Written or the ID case, the Laplacianoperator,.,has the following form(ftomwhrchthe3D- 





(•■1 


...1-2 1 0 0... 






...0 1-2 1 0... 






...0 0 1-2 1... 











(10) 



,00981 Equation ,0 above shows the Laplacian operator for the ID case, because 
U^eLaplacian operator canbe easily shown asamatrix acting onavectorforftelDcase. in 
igLdimensillitiestheLaplacianoperator also c»be expressed asama„x.b. 
Itementsofneighboringpointsin^DorSDspaccc^otberepresented.ana^^^^^^^^ 
«,ematrix will looUiUeasparse matrix, ratherthan^^enear diagonal onemthecaseon^ 
geometry. Forexample, if thepointsinSDspacearenumberedasintheexampefH^^^ 
IS.thentheportionoftheLaplacianoperatoractingonthesepointsandcenteredatthepomt., 

will be as the following: 
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Lx= X2-2X1+X3+ Y2-2Y1+Y3+ Zr2Zi+Z3 
+ X7-2X1+X6+ Y7-2Y1+Y6+ Z7-2Z1+Z6 
+X5-2X1+X4+ Y5-2Y,+Y4+ Z5-2Z1+Z4 



1— rrrrritr^^^^ 

observed *an mat compared with a>eCGmemcd applied in 2D case. 

100991 The regularized analog of Equation 3 has the fonn: 

min(p|P,.-/,llnl|i-ll] (11) 
..ere^lol^egulari^tionparameter. ..e— dinganalos 

new, regularized, matrix A of the form 

ofcircuitpnntedboards. The example ... ^, primitives are combined into a compact 

c Ai «;n9 a cvlinder 503, and a box 504. The pnmiiwc:, oi 
spheres 501 and 502, a cyimaer d , , sn2 and cvhnder 503. But, as 

• wr^. w 504 nartially penetrates into spheres 501 and 502 and cyimae 
group m which box 504 partially P „ea is taken from the corresponding bodies, 

describedbelow,across-sectionmtheoverlappedareaist^^^^^^^^^^^ 

notfromtheboxinthisexample. The scene's volume consists of 32 
example. 
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. 1 t;.n. <;i7es and X-ray cross-sections are shown in the dialog 
101021 The exact locations, sizes, ana ^ , . ^ 

ball 501 is assigned a rank of 3, and ball « as 

definesthecross-section for the volume Where two or^orepnnnj^^^—^^^^^^^^^^ 
cross sectionistakenfromthehighestrarfcedprinuuve. For example, mstde the vo 
r5oTLectsbox504.theorosssectio„ofhal,50Ustakenastheerosss.o„ofthe 

vre:eea„sesuehhall501hasthe(s„perior,r^3,wHilehox504h.r^l. 

1,103, ^esceneSOOisiU— dbypoint-liUesot^sSOSarr^mgedatthe^^^^ 
„,..J,.i.so.^.OSinthisex.ple,..heX-.s^^^^^^^^ 
heamsgoinsinthedireetion.m.e.m^^^^^^ 

-.a^gde^ee......^^^^^^^^^ 

illumination angles varymgbe.ween00a^2,^^;5„,,^„,^^fl,eproblen, handling 
presents .he cases withballgnd arrays (BGA^.=« >»x^ '^^^^^^^^^ 

anobiecthavingparalleltopandbottomsurf^es. F^^^P^,^^^^^^^ ^^^^^^ 

„f.esoeneintheX.dVdiree«o..sen..y.^^^^^^ 
positionforthetopandbottomsurfacesofthebox. Inour , 

X.Yspen,bu, it still hasasizebig^ough to challenge the algona™. 
used to improve the visibility of details in the picture. 
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101051 The goal of 3D tomography is to reoomtruot 3D image data (e.g., 3D cross- 
sections, such as shown in tile 700 of FIGUKE 7) .on, a set of 2D X-tay images cap^^^^ 
regist^tionplane 506. inthis test geometry,.he2DX.ray images .e a. showmF ™ 

, 1- «,ot «f9n Y-rav imaees obtained from 16 ditterent 
and 8B FIGURES 8A and 8B each shows a set of 2D X-ray images 
:ltionsofX-ray.eams0.di«ere„t.Wcsofvie«"or"fie.<.0f.ew",hy.^^^^^ 

Plane 506 FIGURE 8A is shown in the linear scale, and FIGURE 8B .s shown after the log- 
:!::!lforn.tionHghlightingtheproiectionofI«x504.whichispoor.yvisih.em^^ 

scaled images of FIGURE 8A. 

10106, taapreferredembodimen,.heCGalgorithmisusedforproce^ing*e 
eapt»ed2DX-rayimagesforreconstructing3Dimagestheref^m.AsdescrihedaW^^^^^ 

a. orithmisar-iterativealgorithn. ThegraphofnGURB9showsthepro^ss^^^ 
algorithmaftereachstepofiteration-ToexpIaintheourvesmthegraphofFGI^S^ 
excurseintoancxamp.eofacma,implem».ationoftheCGmeth«l.seonstdered.The 
^pseudo-codcdescribesanex^pleimplementationoftheCGalgortthmmaccordance 

with an embodiment of the present invention: 

Compute ro=do =b-A*xo for an initial guess xo 
for j=0,l,..,until convergence 
aj=(rj,rj)/(Adj ,dj) 

X3+i=Xj + ttj dj 

rj+i=rj-aj Adj 
Pj=(rj,rj+i)/(rj,rj) 
dj+i=r Pj dj 
endfor 

whereOisthedotproduct. Two auxiliary vectors, r, d. have the same size as vector X. Due 
to the fact that the exact solution X. (in otherwords: the cross-sections, or vo,ume)is not 
neshouldtesttheco„verge„ceindirect.y.Usua.Iythetestish.edonthcva..of*er^^^ 

vector r. So. curve 901 on graph 900 of FIGURE 9 is e,uM to (r^). Curve 902 shows the 
average square differencebetween actual 2DX-ray images. f,and.he2DX.ray images 
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p^oduoedby me resto«dx:'Z(i',-/.)'-C-e 903 shows «.e average squared difference 

be«.ee„d.e.sh..deof*e:;ua«o„and*eresU>redX,i.e.(Ax.b)^.A„dfinal.y.^ 
repr^en^me average ditfere„oebe.wee„exac.solution,WwWeMsknow„ for memode^ 

and the current x ftom the CG method, i.e. (x«a«- • 

[01,7, Lisa known problem of the CG method when one attempts to use it for the 
2D tomo^aphic restoration thatthemethoddivergesinxspace. as recognizedbyElenaU>h 

P^lm^andFabianaZamain'-TheConiugateGradientRegulariz^^^ 
«Z,Problems^Depar.mentofMathema.ics,Universit.ofB„.ogn^ 

-^««-vo..t02,pagesS7^^^ 
tbeenincor^«tedhereinby.eference.h.viewofthe above, wehaverecogn.edt^^^ 

divergencemamfests itself when the CGmethod is usedin3D.omography.as wen. 

10,08, Becausethemethodconvergesfairlywellonthemodelgeometriesinthe 
spa«ofX-rayimages(curve902),.hesimultaneousdivergencei«Xspaceisiusta 

Lfestationofthefactthatthelinearsystemisill-posed. ^ 
n„mberofX.rayimagesinvolved,the.isstillnoguaranteeofconvergencemxspac e. Tl.».a 

IlonintheCGmethodbecomes: couldonegettheconvergenceof cross-sect— ga 

T \ set of 2D X-ray images obtained ftom different directions (i.e., obtamed 6om 

r:rgir«B:.e:enifsuchasete.sM.couldre,n^ 

toe to be practically useful for the industrial mspectionenvuonment(e.g., for mspectmgP 

Td/orothLevices, That is, too many iterationsoftheCGmethodmayberecm^ for 

:".heconvergencetoenabletheCGme.hodtobepracticallyuseMmh.gh-spee^^^^^^ 

:! ere onstructio„enviromnents,suchasmanyind„st.ialinspec.ionenv— 

Jcalapplicationsthatneedhigh-speedprocessing. The major mamfestat,on o^^e 
:Lgencelthefollowmg..as.henumberofiterationsprogresses,theCGalgor.tlm^^^^ 
migrate into areas withnegative values in vector X,whichhasnophysioal sense mdeed.(n>e 

problem is also observed in the ART family of iterative methods). 

[0,09, Theexample^aphofFlGURE ,0 demonstrates the observed migration of 
«,esolu Jintonon-physieaKnegativOeross-sections. Horizontal Une 1001 ongraph lOOOof 
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FIGURE 10 shows the zero level. The values below such line 1001 are non-physical. 
Unfortunately, it appears that positiveness cannot be imposed on the vector x by a linear 
operator. 

[0110] Choosing the regularization parameter, y, sets the balance between 

minimizing the residual norm 'f|| P,x - f, \\ and mimmizing the "roughness" H Lx \\ {See 

1=1 

Equation 1 1 above). A convenient graphical tool for setting y is the "L-curve". When one plots 
log Ell - fi II ^^^"^^ II ^ II a characteristic L-shaped curve with a (often) distinct 
comer separating vertical and horizontal parts of the curve is obtained. The rationale for using 
the L-curve is that regularization is a trade-off between the data misfit and the accuracy of 
solution. In the vertical part of the curve, the solution norm 'f |1 P,x - || is a very sensitive 

function of the regularization parameter because the solution is undergoing large changes with y 
in an attempt to fit the data better. At these low levels of filtering by the smoothing matrix, L, 
there are still components in the solution which come fi-om dividing by a small singular value in 
the matrix X^/^/ (^^^ Equation 12 above), which corresponds to having inadequate 
regularization. 

[0111] On the horizontal part, the solutions are not changing by very much as y is 
changed. However, the data misfit increases sharply with more filtering. Along this portion, the 
preference for the more highly filtered solutions increases only slightly at the cost of a rapidly 
rising data misfit. FIGURE 1 1 shows graph 1100 that shows an example L-curve that is 
associated with the tomographic restoration in the CG method of an embodiment of the present 
invention. More specifically, FIGURE 11 shows an L-curve for the CG-method calculated for 5 
iterations for the test geometry (of FIGURE 5), and covering the regularization parameter, y, 
between 0.01 and 1000. For each regularization parameter, along the Y-axis of graph 1100 the 
corresponding log of the regularization term is plotted, and along the X-axis of graph 1 100 
different norms of convergence are plotted. 



25237482.1 



36 



Docket No.: 10020666-1 



10U21 There are four graphs plotted in graph 1 100 which show the L-curve for 
different values. Curve 1101 represents the residual, r (as with curve 901 of FIGURE 9); curve 
1 102 shows the difference between actual 2D X-ray images,/ and the 2D X-ray images 
produced bythe restored x:'|(P,x-/,)Mas with curve 902ofFIGURE9);curvell03 shows 

the difference between the right side of the equation and what comes from the restored x, i.e. 
(Ax-b)^ (as with curve 903 of FIGURE 9); and finally, curve 1 104 shows the difference 
between the exact solution, x„..,. which is known for the model, and the solutionxfrom the CG 

method, i.e. (x„,«- X)^ (as with curve 904 of FIGURE 9). If curves 1101. 1102, and 1103 
behave in accordance with the classical L-curve. curvell04,is kind of symmetricallyreflected 
due to the observed divergenceoftheCG method in the space of restored cross-secfons. 

101131 As experiments show, the number of "waves" on the L-curve 1 104 is equal 
to the number of iterations for the CG method: 5 iterations in the example plotted case of 
FIGURE 11 with L-curve 1104 thus having 5 waves (labeled A, B,C,D, and E, respectively m 
FIGURE 1 1). The best visual quality of restoration of the cross-sections corresponds to y-l.O, 
which in turn, corresponds to the first knee calculated torn the left to the right side of graph 

1,00 for curves 1101. 1102, and 1103, or for the first knee calculated from the right side of 
graph 1100 for L-curve 1104. Graph 1100 covers the span of y between 0.01 and 1000. 

10114] As described fiMher below, the quality of the resulting 3D images using the 
above-described CGmethod is very good, generally beingmuch better quaUty than achieved 

through tomosynthesis. While the observed divergence is a problem for a big number of 
iterations, an advantage of the CG method is its superior 3D resolution in cross-secfion m>ages 
achieved just afier very few iterations (~5) compared with the ART technique and tomosynthests 
(e g traditional techniques not using exact solutions of equations). So, by properly choosmg the 
nmnl,erofiterationstobeperfonned,thedivergenceproblemcanbealleviated. And,alarge 
numberofiterations is typicallynot desired in any case, because are computation time is dtrectly 

proportional to the number of iterations. 

[0115] FIGURES 12A-12B and 13A-13B provide a visual comparison of the 
results obtained through the CGmethod(after5iterations)and the results obtained via 
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„*esis«=*a. doe. not use an exact solutionofc^uaaons forme example so«^ 

of FIGURES 5A-5D. More .peoifically, FIGURE 12A show, a tile 1200 .hat compnses 32 

. , ■ Trm re F 12R shows tile 1201 that compnses 
sUces resulting ftom a tomosynthesis techmque. FIGURE 12B shows u 
the32slicesofFIGURE 12A with logteightnessnonnalizationhavingbeenapphed thereto. 

FIGURE 13A shows a tile 1300 that comprises 32 slices .suiting ftom an example CG 
processingmethodofanembodimentofthcpresent invention. andFIGURE13Bshows^e 

:Lcomprisesthe32slicesofFIGURB13Awithloghri^^^^ 

Llicd «teJ. Tiles 1200, 1201. 1300. and 1301 each represent a 3D o^^«te sc.e J 

omGURES5A-5Dby slices talcenperpendicularto the Z-axis at different hetghts alongtheZ- 

a^rThe32slicesare„„merated0eft.bot.om)ineachtile,»darecountedftommehottom^ 

rvolume(i.e..s.icelcor.spondstoabottomcrosa-sectio„belowslice2a,ongtheZ-axts.and 

SO on). 

,01161 As can be seen in this example, the 3D images resulting ftom the CG 
processingteohniquearemuchsharperanddonothave.hesmoothspatialvariations.nte 

:3cLUons typical forthetomosynthes.stechni,ue.Theresol„tionin*^^^ 
^esultingftomteCGprocessing technique is alsobetterti,aninte.omosyn.hes,s techmque. 

,01171 Turning to FIGURE 14. an example hi^-level block diagram is shown tor 
asystemimplementingaJDimagere^nstruction process of an embodiment oftep^^^ 
invention, ht the cxampleofFIGUREM. inspection systemUOOcompnsesaradograpluc 
Lgi„gsysteml402,whichmay.for example. beaco„e-beamtomo^ph.c.mag,ngsysem 
ZasdescribedabovewititFIGUREl. Radiographic imaging system 1402 ,s preferably a 
Cl-:imaginssystemtitatcapn^20radiographicimagesofano.ectl^^^^ 
JaniyReconstructionprocessorUOSiscommunicativelycoupledtoradto^^^^^^^^ 
systeml4^suchti.at it canreceive 2D radiograpMc image da.a(or>xels") ftom .magerl«^^ 
Lns«uctionprocessorl403 is ako communicatively coupled to data storage ^r -memory) 
I 04whichis:ed.ostorerecons.ructed3D(e.g..cross.sectional)imagedata. D.^^^^^^^ 

1404maycompriseanysuitabledatastoragedevicenowlmownorlaterd.scove^for^^^^^ 
_nstructed3Dimage data. includingas examples random access mem^(RA^*^^ves. 

floppy disks, optical discs(e.g..Compao.Discs(CDs)andDigital Video D.SCS (DVDS)), and 
!rdatast;rrgedevices.DaUstoragel404maybeintegratedwiti,reconstruotionprocessor 

1403 (e.g., on-chip memory) or it may be external ttiereto. 
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10118] As described above, in accordance with embodiments of the present 
invention :«=onstr«ction processor 1403 may implement an iterative processing technique, such 
as a CG method, for processing the 2D image data to reconstruct 3D image data therefrom. 
Reconstructionprocessor 1403 may be implemented as the reconstruction processor disclosed >n 

co-pendingU.S. Paten. Application SerialNumber [Attorney Docket No. 10020619-1] 

titled "A SYSTEM AND METHOD FOR PARALLEL IMAGE RECONSTRUCTION 
OF MULTIPLE DEPTH LAYERS OF AN OBJECT UNDER INSPECTION FROM 
RADIOGRAPHIC IMAGES", the disclosure of which is hereby incorporated herem by 
reference. Accordingly, in certain implementations, reconstruction processor 1403 maybe 
implemented to process multiple depth layers of object 10 in parallel, such as disclosed m U.S. 

Patent Application Serial Number [AtU>mey Docket No. 10020619-1] titled "A 

SYSTEM AND METHOD FOR PARALLEL IMAGE RECONSTRUCTION 
OF MULTIPLE DEPTH LAYERS OF AN OBJECT UNDER INSPECTION FROM 
RADIOGRAPHIC IMAGES", whereto an iterative processing technique, such as a CG 
algorithm maybe used for performing suchprocessingtoaccordance with embodimentsofthe 

present invention. Further, in certain implementations, reconstruction processor 1403 may 
acmallycompriseapluralityofprocessors that executemparallel for reconstructingSDmrages 

of object 10 from 2D radiographic m,ages received from imager 1402, whe^n such pluralrty of 
processors may each utiUze an iteradve processing technique, such asaCG algorithm for such 

processtog m accordance with embodiments of the presem invention. 

[0119] Embodiments of the present mvention may be implemented as a 
combmation of hardware and firmware (or circuitry) that receives as input 2D radiographic (e.g.. 
X-raytr«sive)image data for at leastapordonofan object and recons,ructs3D(e.g., cross- 
sectional) image data corresponding therefrom using an iterative processing technique, such asa 
CG algorithm. A preferred embodiment comprises a hardware architecture that includes a 
processor unit such as a Field Programmable Gate Array (FPGA) or an Application Specific 
integrated Circuit (ASIC). Thus, reconstruction processor 1403 ofFIGURE 14 is preferably 
implementedaspartofalargerprocessorumtsuchasanFPGAorASIC. Variousother 
configurations now known or later developed for implementingareconstruction processor for 

processtog 2D mrage data to reconstruct 3D image data therefrom to accordance wrth 
embodiments of the present invention may be used m alternative embodiments. 
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[0120] In view of the above, in operation of a 3D tomographic system 
implementing the CG method, one may implement a CG method for iteratively processing 2D 
image data for reconstructing 3D image data therefrom. More specifically, one may select a 
regularization parameter (which may be determined from the L-curve in the manner described 
above), select an optimal number of iterations to be performed for the particular imaging 
geometry and/or object being imaged by such system (wherein such optimal number of iterations 
may be selected from a lookup table of previously determined iterations that are optimal for a 
given type of imaging geometry and/or object being imaged), and then use the CG method for 
iteratively processing 2D images received from a radiographic imaging device to construct 3D 
images. 

[0121] Although the present invention and its advantages have been described in 
detail, it should be understood that various changes, substitutions and alterations can be made 
herein without departing from the spirit and scope of the invention as defined by the appended 
claims. Moreover, the scope of the present application is not intended to be limited to the 
particular embodiments of the process, machine, manufacture, composition of matter, means, 
methods and steps described in the specification. As one of ordinary skill in the art will readily 
appreciate from the disclosure of the present invention, processes, machines, manufacture, 
compositions of matter, means, methods, or steps, presently existing or later to be developed that 
perform substantially the same fimction or achieve substantially the same result as the 
corresponding embodiments described herein may be utilized according to the present invention. 
Accordingly, the appended claims are intended to include within their scope such processes, 
machines, manufacture, compositions of matter, means, methods, or steps. 



25237482.1 



40 



